Constraints on Modified Gravity from ACT and SPT 
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The Atacama Cosmology Telescope (ACT) and the South Pole Telescope (SPT) have recently 
provided new and precise measurements of the Cosmic Microwave Background anisotropy damping 
tail. This region of the CMB angular spectra, thanks to the angular distortions produced by grav- 
itational lensing, can probe the growth of matter perturbations and provide a new test for general 
relativity. Here we make use of the ACT and SPT power spectrum measurements (combined with 
the recent WMAP9 data) to constrain f(R) gravity theories. Adopting a parametrized approach, 
we obtain an upper limit on the lengthscale of the theory of Bo < 0.86 at 95% c.l. from ACT, while 
we get a much stronger limit from SPT with Bo < 0.14 at 95% c.l.. 

PACS numbers: 98.80.Es, 98.80.Jk, 95.30.Sf 



I. INTRODUCTION 

The major goal of modern cosmology is to understand 
the source of cosmic acceleration. One of the possible 
solutions to this puzzling phenomenon is to modify gen- 
eral relativity on very large scales in order to allow an 
accelerating phase in matter-only universes. Examples 
of such "Modified Gravity" (hereafter MG) models are 
f(R) theories [1-5] and in the recent years several au- 
thors have searched for modified gravity and departures 
from general relativity in cosmological data |6l-[2l| 

The recent precise measurements of the Cosmic Mi- 
crowave Background damping tail from the Atacama 
Cosmology Telescope (ACT) [H| and the South Pole 
Telescope (SPT) [23[ are offering a new opportunity to 
further test MG theories. 

The shape of the damping tail of the CMB anisotropics 
depends strongly from the effect of lensing caused by the 
intervening matter densities along the line of sight of the 
CMB photons. CMB lensing therefore probes the growth 
of perturbations up to redshift z ~ 6. Since the ampli- 
tude and the evolution of matter perturbations can be 
drastically altered in MG theories a precise detection of 
the CMB lensing can place strong constraints on these 
deviations and possibly identify them (see e.g. [Hj]). 

However the ACT and SPT experiments are reporting 
quite different constraints on the amount of CMB lens- 
ing(see the discussion in 24]). Parametrizing the lensing 
amplitude by an effective amplitude Al, that is Al = 1 
in case of the standard exp ected signal and Al = in 
case of no lensing (see [26| for a definition), the ACT 
data provide the constraint Al — 1.7 ± 0.38 at 68% c.l. 
([22|). therefore indicating a larger amplitude, while the 
SPT is more consistent with the standard expectations 
with A L = 0.85 ± 0.15, again at 68% c.l. ((H). 

The Al parameter is clearly an effective parameter 
and can be used just to indicate possible deviations from 
the expectations of the standard scenario. It is there- 
fore timely, as we do in this paper, to analyse the results 
from ACT and SPT in the context of more physically 
consistent scenarios, as MG theories. 

Here we adopt the parametrized modified gravity sce- 



nario presented in [4f| , restricting our analysis to the case 
of f{R) theories. In this model, the background expan- 
sion is identical to the one produced by a cosmological 
constant, while the evolution of perturbation is altered 
and depends on a single parameter Bq that represents 
the length-scale of the theory [27j . 

Since the ACT and SPT datasets are providing results 
that are significantly different, we take a conservative 
approach to analyse each dataset and discuss the corre- 
sponding results separately. The ACT and SPT datasets 
are combined with the recent data release from nine of 
observations from the Wilkinson Microwave Anisotropy 
Probe (WMAP9) [H, 

The paper is organized as follows. In Section [TTJ we 
present the modified gravity model considered for our 
analysis, in Sec. IHII we describe the analysis method and 
in Sec. II VI we present our results. We conclude in Section 
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II. PARAMETRIZED f(R) GRAVITY 

The f(R) theories are currently one of the most pop- 
ular class of MG models. These models generalize the 
Einstein-Hilbert action replacing the Ricci scalar with a 
function of R itself. The generic modified action is 



S = 



■W) 



(1) 



where k 2 = 8irG (c=l) and C m is the matter la- 
grangian density. 

Focusing on this particular MG category is interesting 
for two main reasons. Firstly, their modified Lagrangian 
is quite simple and generic, since the modified dynamic at 
every scale is recovered using only the first order invari- 
ant. Secondly, some models belonging to this class have 
been shown to satisfy both cosmological viability condi- 
tions and local tests of gravity, thanks to the chameleon 
mechanism [37H39l | . 

In order to reproduce the effects of f(R) gravity in the 
evolution of matter perturbations here we adopt a generic 
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MG parametrized form, proposed in [40(, specializing it 
to the f(R) case. In this parametrization the background 
is fixed to that of ACDM and the modifications in the 
linearized Einstein equation are encoded in two scale- and 
time-dependent parametric function p,(k, a) and j(k, a) 



fc 2 * = -fi(k, a)4TrGa 2 {pA + 3(p + P)a} 



(2) 



k 2 [<S> - 7 (jfe, a)*] = a)l2irGa 2 (p + P)a (3) 

where $ and $ are the two scalar metric potentials in 
the Newtonian gauge, cr is the anisotropic stress that 
vanishes for baryons and CDM but not for relativistic 
species, 5 = Sp/p is the density contrast and pA is the 
comoving density perturbation, defined as 



Ha , „. 
pA = pS + 3—(p + P)v 

k 



(4) 



where v is the velocity field. 

It has been shown in [4l[ that we can recover the 
f(R) theories choosing the following parametric form for 
p(k, a) and j(k, a) 



2.4 km s 1 Mpc 1 , consistently with the measurements of 
the HST [H. 

We include information from measurements of bary- 
onic acoustic oscillations (BAO) from galaxy surveys, 
combining four datasets: 6dFGRS from [3(J, SDSS-DR7 
from SDSS-DR9 from [H and WiggleZ from [35j]. 
We refer to this dataset as BAO. 

We sample a seven-dimensional set of cosmological pa- 
rameters, adopting flat priors on them: the Bq modi- 
fied gravity parameter, the baryon and cold dark matter 
densities Q\,h 2 and fl c h 2 , the ratio of the sound horizon 
to the angular diameter distance at decoupling 0, the 
optical depth to reionization r, the scalar spectral in- 
dex n s , the overall normalization of the spectrum A s at 
k = 0.002 Mpc" 1 . 

Given the tension between the ACT and SPT exper- 
iment in the lensing ampllitude, we also consider varia- 
tions in the lensing amplitude parameter Al as defined 
in 26] . Finally, the amount of helium abundance in the 
universe Y p is fixed by assuming Big Bang Nucleosynthe- 
sis in the standard case of three neutrino families. 



IV. RESULTS 



p(k,a) 



1 + f A 2 k 2 a s 



1 



X 2 k 2 a s 



7(/c,a) 



1 + §A 2 k 2 a s 
1 + fA 2 k 2 a s 



(5) 



Viable f(R) models must have s ~ 4 in order to closely 
mimic ACDM expansion (27j, that is the case we are 
interested in. Indeed the only free parameter we consider 
in our analysis is the characteristic lengthscale Ai. It is 
usual expressed in literature in term of the dimensionless 
parameter Bq as follows: 



Bq = 



2H 2 \\ 



(6) 



i.e. it gives the lenghtscale in units of the horizon scale. 



III. DATA ANALYSIS METHOD 

Our theoretical models are computed with the public 
available code MGCAMB (ibj v .2 while the analysis is 
based on a modified version of CosmoMC [2!| a Monte 
Carlo Markov Chain code. 

We consider the following set of recent CMB data 
(publically available on the corresponding web pages): 
WMAP9 [H], SPT [13, ACT [H including measure- 
ments up to a maximum multipole number of i max = 
3750. 

For the ACT experiment we use the "lite" version of 
the likelihood (3(| that has been tested to be correct also 
in the case of the extension respect to A-CDM models. 

We also consider a gaussian prior on the Hub- 
ble constant (hereafter HST prior) Hq = 73.8 ± 



Our main results are reported in Table I. Since the 
ACT and SPT datasets are reporting significantly dif- 
ferent constraints on Bq we consider these two datasets 
separately. 

As we can see, both ACT and SPT are not providing 
any evidence for MG. However, the SPT dataset gives 
significantly stronger constraints on Bq (Bq < 0.14 at 
95% c.l.) respect to those derived by ACT (B < 0.90 at 
95% c.l.). The difference appears as even more striking 
in Figure 1 (Left Panel), where we report the two poste- 
riors on Bq coming from the two experiments: while SPT 
strongly constrain Bq, the posterior from ACT shows a 
bimodal distribution, suggesting an higher compatibility 
with modified gravity models. The reason of this differ- 
ence is mainly due to the differenr lensing signal present 
in the ACT e SPT TT spectra (see H): since f(R) MG 
models increase the lensing signal they are more consis- 
tent with the larger amplitude of ACT than with the 
smaller amplitude of SPT. The best fit value for ACT 
is indeed Bq ~ 0.78 even if this dataset still does not 
provide any compelling evidence for MG. 

The inclusion of the HST prior and of the BAO dataset 
improves the constraints (Bq < 0.12 at 95% c.l. from 
SPT and B < 0.86 at 95% c.l. from ACT), however 
not in a significant way, clearly showing that most of the 
constraining power is coming from the CMB spectrum 
distortions introduced by gravitational lensing. 

It is interesting to consider the impact of MG on the 
standard cosmological parameters. As we see from the 
Table, and as already showed in L 40j, there is little cor- 
relation between Bq and the, standard, six cosmological 
parameters. We found that the largest correlations are 
with scalar spectral index ns and amplitude A$. How- 
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FIG. 1. Posterior distribution functions for the -Bo parameter in the case of Al = 1 (Left panel) and Al free (Right Panel). 
The different lensing amplitude measured by ACT makes MG model more consistent with the data and a bimodal posterior 
distribution is present (Left Panel). When variations in Al are considered the SPT bound is weaker, while the ACT dataset 
is more consistent with GR. The bimodal distribution still present in the right panel is due to the low WMAP anisotropy at 
large angular scales. 
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FIG. 2. Effect of Bo on the CMB angular temperature power spectrum. We plot the differences respect to the standard A-CDM 
model. On the left panel we see that the anisotropy at low multipoles decreases as Bo is increased. A larger Bo is therefore 
more compatible with the low WMAP quadrupole. On the right panel we see that the anisotropies on large angular scales are 
increased respect to A-CDM. The effect is due to an increase in the CMB lensing amplitude. 



ever, these correlations changes in function of Bq. When 
Bq << 1, larger Bq is in more agreement with smaller 
ns and larger As- When Bo ~ 1, larger Bq is in more 
agreement with larger ns and smaller As- 

In order to further test the importance of the lensing 
signal in constraining modified gravity models we have 
performed an analysis by letting variations in the lensing 
amplitude A L . The results are reported in Table 2. As 
we can see, the effect of marginalizing over the lensing 
amplitude is clearly to make weaker the SPT constraint 
and to leave as unaffected the ACT constraint. However 
the bimodal distribution present in the ACT case is now 
suppressed as we can see from Figure 1 (Right Panel) 
where we plot the posterior distribution functions for Bq. 
Moreover, as we can see from the results in Table 2, the 



ACT lensing signal is consistent with A L = 1 in the case 
of MG models. We can therefore conclude that while 
the ACT data does not show any evidence for MG, the 
lensing signal is in better agreement with the A L = 1 
case in the framework of MG. 



It is interesting to note that, while now suppressed 
from the previous case, the bimodal distribution is still 
present when A L varies and it is also now evident in 
the SPT dataset. The reason is that MG f(R) gravity 
models produce also a lower quadrupole and lower £ tem- 
perature anisotropy in agreement with the WMAP data 
(see Figure 2, left panel). 
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Parameters 


SPT+WMAP9 


ACT±WMAP9 


SPT±WMAP9±HST±BA0 


ACT±WMAP9±HST±BA0 


tt c h' 2 

woe 

T 

n 3 
B 

B [km/s/Mpc] 
log(10 10 A s ) 

O m 

Age/Gyr 
n sz 

^3000 

n CL 

-^3000 

n PS 

-^3000 

A sz 


0.02224 ± 0.00034 
0.1091 ±0.0036 
1.0428 ± 0.0010 
0.0827 ± 0.013 
0.9676 ± 0.0093 

< 0.14 (95% c.l.) 

72.2 ± 1.7 
3.060 ± 0.027 
0.747 ±0.018 
0.253 ±0.018 
13.689 ± 0.066 

4.2 ± 2.1 
4.8 ±2.0 

20.3 ± 2.4 


0.02281 ± 0.00044 
0.1142 ± 0.0044 
1.0402 ± 0.0020 
0.091 ±0.014 
0.973 ±0.012 
< 0.90 (95% c.l.) 
70.2 ± 2.1 
3.174 ±0.045 
0.721 ±0.025 
0.279 ± 0.025 
13.71 ± 0.10 

0.94 ±0.57 


0.02235 ± 0.00033 
0.1119 ± 0.0023 
1.04257 ± 0.00098 

0.080 ±0.012 
0.9633 ± 0.0078 
< 0.12 (95% c.l.) 
70.9 ± 1.0 
3.066 ± 0.025 
0.733 ±0.012 
0.267 ±0.0012 
13.724 ± 0.053 
4.0 ± 2.1 
4.8 ±2.0 
20.5 ±2.3 


0.02279 ± 0.00040 
0.1149 ±0.0028 
1.0403 ±0.0019 
0.090 ± 0.013 
0.9724 ± 0.0096 
< 0.86 (95% c.l.) 
70.0 ± 1.3 
3.185 ±0.035 
0.718 ±0.015 
0.282 ± 0.015 
13.714 ± 0.084 

0.91 ±0.56 


Xmin/2 


3808.25 


3799.09 


3811.41 


3800.89 



TABLE I. Constraints on the MG parameter Bo and the standard cosmological parameters described in the text from ACT 
and SPT combined with WMAP9, HST prior and BAO. We report constraints at 68% confidence level (bounds on Bo are at 
95% c.l). 



Parameters 


SPT±WMAP9 


ACT±WMAP9 


SPT±WMAP9±HST±BAO 


ACT±WMAP9±HST±BAO 


n c h' 2 

9 

T 

A L 

n 3 

Bo 

B [km/s/Mpc] 
log(10 10 A s ) 

Age/Gyr 
n sz 

-^3000 
n CL 

n PS 

A S z 


0.02204 ± 0.00033 
0.1154 ±0.0027 
1.04200 ± 0.00097 
0.083 ±0.013 
0.60 ±0.10 
0.9561 ± 0.0084 
< 0.73 (95% c.l.) 
69.1 ± 1.2 
3.083 ±0.026 
0.712 ±0.015 
0.288 ±0.015 
13.790 ± 0.055 

5.1 ±2.3 

5.2 ±2.1 
20.0 ±2.4 


0.02294 ± 0.00048 
0.1132 ± 0.0043 
1.0406 ± 0.0019 
0.090 ±0.014 

1.25 ±0.29 
0.971 ±0.011 
< 0.91 (95% c.l.) 
70.8 ± 2.1 
3.177 ±0.040 
0.727 ± 0.024 
0.273 ± 0.024 
13.68 ±0.10 

1.9 ±1.3 


0.02215 ±0.00032 
0.1137 ± 0.0025 
1.04233 ± 0.00096 
0.085 ±0.013 
0.62 ±0.11 
0.9598 ± 0.0081 
< 0.77 (95% c.l.) 

70.0 ± 1.1 
3.082 ± 0.026 
0.722 ±0.013 
0.278 ± 0.0013 
13.760 ± 0.053 

5.0 ±2.3 
5.2 ±2.1 

20.1 ± 2.4 


0.02289 ± 0.00039 
0.1144 ±0.0027 
1.0403 ±0.0018 
0.089 ± 0.013 
1.19 ±0.26 
0.9699 ± 0.0097 
< 0.85 (95% c.l.) 
70.2 ± 1.2 
3.183 ±0.034 
0.721 ±0.014 
0.279 ±0.014 
13.702 ± 0.079 

1.7 ± 1.2 


Xmin/2 


3807.34 


3798.95 


3808.98 


3800.64 



TABLE II. Constraints on the f(R) parameter Bo, the lensing amplitude Al and the standard cosmological parameters 
described in the text from ACT and SPT combined with WMAP9, HST prior and BAO. We report constraints at 68% 
confidence level (bounds on Bo are at 95% c.l). 



V. CONCLUSIONS 

In this brief paper we have presented new constraints 
on f(R) MG models from the new recent measurements 
of the CMB damping tail provided by the ACT and SPT 
experiments. We have found that both experiments show 
no evidence for deviations from GR. However, while the 
SPT data significantly improves the previous constraints 
obtained from similar analysis, the ACT data gives much 
weaker constraints and shows a bimodal posterior dis- 
tribution for Bq. We have attributed this different be- 
haviour to the different amplitude of the lensing signal 



detected by those experiments and showed that when 
the lensing amplitude A L is let to vary both datasets 
provide similar constraints. When A L is varied, we have 
found that the ACT data does not show any indication 
for Al > 1 in the framework of MG models. Moreover, a 
bimodal distribution for Bq is present in both ACT and 
SPT datasets when we marginalize over A L . This is due 
to the large angular scale regime of the measured CMB 
spectrum, that prefers a low quadrupole and a bluer spec- 
tral index (see e.g. [3^|). 

Presenting combined results from ACT and SPT spec- 
tra as in [44J needs to be done with great care: while 



5 



compatible in between two standard deviations in the 
case of a standard six parameter analysis, the two exper- 
iments could show very different constraints in extended 
theoretical frameworks, especially when the lensing sig- 
nal plays a significant constraining role. 

It is useful to compare the SPT results with previous 
limits on Bq present in the literature. The constrain we 
obtain from WMAP9+SPT+H0+BAO, in the case A L = 
1, is much tighter than the constraint of Bq < 0.42 (95% 
C.L.) obtained from a combined analysis of cosmological 
data and Integrated Sachs Wolfe data [42[ and from the 
similar constraints Bq < 0.4 (95% C.L.) obtained in [|(| 
and B < 0.42 (95% C.L.) from Q, where the datasets 
considered are slightly different between papers. In [HI 
they also found a very tight constrain combining with 
cluster abundance data (Bo < 0.001 95% C.L.), however 
this constraint is obtained in the non linear perturbation 
regime where the simple treatment of f(R) models we 
adopt here may not be sufficient. 

While the SPT provides a much better constraint, 
one should however also consider it with great caution, 
given the tension on the lensing amplitude with the ACT 



dataset. 

Finally, the ACT collaboration has provided a determi- 
nation of the lensing amplitude also from the four points 
CMB correlation function (see [22]). This amplitude is 
perfectly consistent with the standard case, however we 
prefer here to do not include this dataset for the following 
conservative reasons: a) we prefer to compare the ACT 
and SPT datasets at the same power spectrum level; b) 
the ACT constraint from higher correlations comes from 
about 50% of the data used in the estimation for the 
power spectrum (the ACT-E dataset). 

The current measurements of CMB lensing will be dra- 
matically improved by the Planck satellite mission, ex- 
pected to release new data by end on March 2013. 
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